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Abstract 

In this paper we address the issue of the Gribov copies in SU{N),N > 2, Euchdean Yang-Mihs 
theories quantized in the maximal Abehan gauge. A few properties of the Gribov region in this 
gauge are estabhshed. Similarly to the case of SU{2), the Gribov region turns out to be convex, 
bounded along the off-diagonals directions in field space, and unbounded along the diagonal ones. The 
implementation of the restriction to the Gribov region in the functional integral is discussed through 
the introduction of the horizon function, whose construction will be outlined in detail. The influence 
of this restriction on the behavior of the gluon and ghost propagators of the theory is also investigated 
together with a set of dimension two condensates. 



Introduction 



In his seminal work [Ij, Gribov pointed out that the quantization of Yang-Mills theories through the 
Faddeev-Popov method is plagued by the existence of the Gribov copies. Even if one imposes a sub- 
sidiary gauge fixing condition in an attempt to remove the gauge redundancy, there still exist field 
configurations belonging to the same gauge orbit which fulfill such condition, i.e. there will be equivalent 
field configurations, or copies, in the gauge fixed theory. As a consequence, the functional measure in the 
Feynman path integral becomes ill-defined. 
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In order to circumvent this problem, Gribov suggested [Ij that the domain of integration in the Feynman 
path integral should be restricted to a certain region, known as the Gribov region, which was supposed to 
contain inequivalent field configurations only. Nowadays, it is known that the Gribov region is not com- 
pletely free from Gribov copies, i.e. additional copies still exist within the Gribov region [21 [3]. In order 
to avoid the presence of these additional copies, a further restriction to a smaller region, known as the 
fundamental modular region, should be implemented [2., 3j. Though, a procedure of effectively achieving 
the restriction to the fundamental modular region is still beyond our present capabilities. Therefore, in 
the following, we shall limit ourselves to investigate the Gribov region. 

In [l], the Landau gauge was employed in order to illustrate several aspects related to the existence 
of the Gribov copies. In subsequent works [H HJ HI [6] , properties of the Gribov region in this gauge were 
established. In particular, in [7,^, Zwanziger was able to show that the restriction to the Gribov region 
could be achieved by adding to the Faddeev-Popov action an additional nonperturbative term, called 
the horizon function. Remarkably, the resulting action, known as the Gribov-Zwanziger action, can be 
cast in local form and proven to be multiplicatively renormalizable [3 [51 O [TOl [TT] . More recently, in 
\12\ [T3 l [T ^ I15j. dynamical effects related to the condensation of local dimension two operators have been 
accommodated in the Gribov-Zwanziger action in a way compatible with locality and renormalizability, 
giving rise to the so called refined Gribov-Zwanziger model. As far as the behavior of the gluon and 
ghost propagators is concerned, the refined model yields a positivity violating gluon propagator which is 
suppressed in the infrared and which does not vanish at zero momentum, while giving a ghost propagator 
which is no longer enhanced in the infrared, behaving essentially as for ~ 0. Such behavior, referred 
as to the decoupling or massive solution, has also been obtained within the context of the Schwinger- 
Dyson equations [HI [T71 [18l [19] . So far, the behavior of the gluon and ghost propagators in the infrared 
region obtained from the refined model seems to be in agreement with recent lattice data in the Landau 
gauge [20\ \21 \ [22l 123 1 I25j . Though, we have to underline that no unanimous consensus on this matter 
has yet been reached, see e.g. [191 [261 [27] for other possible positions. At present, the issue of the infrared 
behavior of the gluon and ghost propagators in the Landau gauge is object of a rather interesting and 
intensive debate. 

The effects stemming from the existence of the Gribov copies turn out to be relevant in the nonper- 
turbative infrared region, and might play an important role for color confinement in QCD. Moreover, 
this is a central issue towards a correct quantization of Yang-Mills theories. As pointed out in [28], the 
existence of the Gribov copies is in fact a general feature of the gauge fixing procedure. Besides the Lan- 
dau gauge, the issue of the Gribov copies has been addressed in other gauges such as: the Coulomb gauge 
[29] and the linear covariant gauge [30], which has the Landau gauge as a particular case. Concerning 
the maximal Abelian gauge, only the particular case of SU{2) has been discussed so far [HI [32] . 

In the present work we shall consider the case of Euclidean SU{N), N > 2, Yang- Mills theories in 
the maximal Abelian gauge. Let us also remind that this gauge is suitable for the study of the so called 
Abelian dominance hypothesis [33 [ I34 [ [35]. which is one of the main ingredient of the dual superconduc- 
tivity mechanism for color confinement [36[ [371 [38] , according to which Yang- Mills theories in the low 
energy region should be described by an effective Abelian theory in the presence of monopoles. A dual 
Meissner effect arising as a consequence of the condensation of these magnetic charges might give rise 
to the formation of flux tubes which confine the quarks. The maximal Abelian gauge displays also the 
important property of possessing a lattice formulation [39 [ I40 [ [^. while being a renormalizable gauge in 
the continuum [i2l [i3l [M [l5] . 

As already mentioned, several works have been devoted to investigate the effects of the Gribov copies in 
the maximal Abelian gauge in the case of SU{2). A study of the influence of the Gribov copies on the 
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two-point gluon and ghost correlation functions was reported in [32j . These propagators have been also 
investigated within the Schwinger-Dyson framework, see [l6]. In [17j, general properties of the Gribov 
region in this gauge were established. In [38], the horizon function was obtained under the requirements 
of localizability and renormalizability. Finally, the inclusion of a set of local dimension two operators in 
the presence of the horizon function was presented in [49|, giving rise to a model which is analogous to 
the refined Gribov- Zwanziger of the Landau gauge. It is worth to point out that the resulting behavior 
of the tree level gluon and ghost propagators obtained in [49] are in good agreement with the most recent 
lattice numerical simulations [41j done in the case of SU{2). 

In the present work, we attempt to generalize the entire program outlined in the series of papers 
[32l|47l|48l|49] to SU{N) Yang-Mills theories in the maximal Abehan gauge. The general case of SU{N), 
N > 2, is interesting for, at least, two reasons. Firstly, as the symmetry group of QCD is SU{2>), the 
case of = 3 is of particular interest. Secondly, as we shall see, for N > 2, the Faddeev-Popov operator 
of the maximal Abelian gauge possesses an additional term which is absent in the case of SU{2). This 
term will induce modifications on the behavior of the correlation functions. For N = 2, only the diagonal 
components of the gluon propagator are affected by the restriction to the Gribov region. However, for 
A > 2, the off-diagonal components are affected too, through a term proportional to (A — 2). This 
feature might raise interest from the lattice community. A numerical study of the gluon propagator in 
the case of SU(3) is, to our knowledge, still lacking. 

The paper is organized as follows. In section 1, the Faddeev-Popov quantization in the maximal Abelian 
gauge is shortly reviewed. In section 2, we introduce the Gribov region and we establish some of its 
properties. This region turns out to be convex, unbounded in all diagonal directions in field space, and 
bounded along the off-diagonal directions. Section 3 is devoted to the proof of a useful statement which 
ensures that for any field configuration belonging to the Gribov region and located near its bound- 
ary, there is a Gribov copy, close to the boundary of fi, however located outside of fl. This statement, 
originally proven by Gribov in the case of the Landau gauge [1], provides us a support for restricting 
the domain of integration in the functional integral to the region $7. The effective implementation of 
such a restriction will be the main subject of section 4. In section 5, in order to go beyond Gribov's 
quadratic approximation, we introduce the horizon function. Although in the case of SU{N) several 
possible candidates for the horizon function can be written down, only one term turns out to be selected 
by the requirements of being localizable, of reproducing the horizon function already introduced in the 
case of SU{2), and of coinciding with the expression obtained from Gribov's no pole condition [T]. In 
section 6, a local action containing a suitable set of dimension two operators chosen in such a way as to 
preserve the symmetry content of the theory is presented. The resulting gluon and ghost propagators 
will be worked out. Finally, in section 7, we collect our conclusion and discuss some perspectives. 



1 The maximal Abelian gauge 
1.1 Some useful definitions 

In order to introduce the maximal Abelian gauge fixing condition, one makes use of the following decom- 
position for the gauge field A^{x): 

Af, = A^T^ = ApP'' + A^p' , (1) 

where are the (A^ — 1) generators of SU{N), while T" and T* stand, respectively, for the off-diagonal 
and diagonal generators. There are A^(A^ — 1) off-diagonal generators and (A — 1) diagonal ones. The 
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diagonal generators give rise to the Abelian subgroup U{1)^ ^, which is known as the Cartan subgroup. 
Here, we have adopted three sets of indices, namely, 

A,B,C,D,E,... G {!,... ,(iV2-l)}, 
a,b,c,d,e,... G {!,..., iV(iV-l)}, 

i,j,k,l,... G {!,... ,(7V-1)}, (2) 

where the capital Latin indices {A, B,C, D, . . . } are the usual SU{N) group indices. The off-diagonal 
indices are represented by {a, b,c,d,e,...} and run from 1 to N{N — 1), while the diagonal indices are 
given by {i,j,k,l, . . . } running from 1 to (N — 1). The commutation relations between the generators 
T^, T'^, r* are given by 

[T\T^] =0, (3) 

where f^^'~^ denote the structure constants of SU{N). From the last equation of ([3|) we conclude that 
jaij _ pjk _ Q rjij^g remaining nonvanishing structure constants f"'^'^ and are totally antisymmetric 
by the exchange of indices and obey the following identitief]: 

Q jabi jbcj _j_ jabj jbic 

Q jabc jcdi _|_ jadc jcib _j_ jaic jcbd 

Q yabc jcde _j_ jabi jide _j_ j^adc jceb _j_ jadi jieb _j_ jaec jcbd j_ jaei jibd 

The off-diagonal and diagonal components of the gauge field transform under an infinitesimal gauge 
transformation as 

6A^^ = -id.co^ + gr'^Ay). (6) 
Here, is the covariant derivative with respect to the diagonal components, 

and (cj'^jO;*) are the gauge parameters. The Yang-Mills action can also be written in terms of these 
diagonal and off-diagonal components as 

Sym = IJ d'x {f;,f;, + , (8) 

with 

F;, = O'l^Al-Dl^A^ + gr^^^A^Al, 

= d^A'.-d.A'^ + gr^AIAl. (9) 



'^These identities follow from the Jacobi identity: 



jABC jCDE _|_ jADC jCEB _j_ jAEC jCBD _ q 
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1.2 Describing the maximal Abelian gauge condition 



Since we have split the gauge field A^{x) into two components, a diagonal (or Abelian) one represented 
by ^^(x) and an off-diagonal one given by we can now choose to fix the gauge invariance of 

expression ([8]) by imposing a gauge condition on each component of the gauge field. For the off-diagonal 
components ^^(a;) one imposes the condition 

DfA'^ = d,A; - gr^'^A^^Al = , (10) 
while for the diagonal ones one requires 

d,Al = Q. (11) 

Note the nonlinearity in the gauge fixing of the off-diagonal components, eq. (jlOp . This particular choice 
of gauge fixing is supported by the interesting property that it follows from the requirement that the 
auxiliary functional 

F[A]= [ d^^xA^Al, (12) 



is stationary with respect to the local gauge transformations ([6]). The nonlinear condition (jlOp still allows 
for a remaining local U{1)^~^ invariance which is removed by imposing a gauge condition on the diagonal 
components, for which a Landau type condition is usually chosen, as expressed by eq. (fTTI) . Conditions 
(llOP and (jlip are referred as to the maximal Abelian gauge. 



For further use, let us also display here conditions (jlOp and (jlip in momentum space. By perform- 
ing a Fourier transformation, one straightforwardly obtains 

k,Al{k) = igr^ j ^^Al{k-p)A'^{p), 

k^A^ik) = 0. (13) 



1.3 The Faddeev-Popov quantization 

The gauge fixing term, which naturally arises from the Faddeev-Popov quantization method, assumes 
the following form in the case of the maximal Abelian gauge: 



•^MAG 



{(^DfAl + c'd^A^^ 



(fx 



+ib' d^A^^ + d^id^c' + gr'^Alc''^ 



(14) 



In this expression, the auxiliary fields {W, ¥) are the Lagrange multipliers enforcing the maximal Abelian 
gauge fixing conditions (fTOl) . (fTT]l : (c",c*) are the Faddeev-Popov ghost fields, while (c",c*) are the anti- 
ghost fields; Ai""^ stands for the Faddeev-Popov operator, given by 



-Dl'D'^ - gr'^AlDl^ - g^nf^^A'^^A'l ; 



(15) 
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and finally, s is the nilpotent BRST operator acting on the fields as 



(16) 

sc" = , sc* = i¥ , 

sb" = 0, sb' = 0. 

The nonlinearity of the gauge fixing implies that an additional self-interacting ghost term is required in 
order to ensure renormalizability [421 B3l [55] . This additional term is defined as follows 



2 2 \ 

+y r^7'"^*c*c^c'^c* + ^r^^f'^^^^^c'^c^ \ . (17) 

Notice that this term is proportional to a gauge parameter a and thus the original gauge fixing is recovered 
when a is set to zero. However, the limit a — )• has to be taken after the removal of the ultraviolet 
divergences. In fact, some of the terms proportional to a would reappear due to radiative corrections, 
even if = 0, [42l|43l|44l|45]. Moreover, the action 

S = Sym + 'S'mAG + Sa (18) 
is multiplicatively renormalizable to all orders of perturbation theory \A2 \ 1431 I45j. 

To conclude this section, let us make a few remarks on the partition function of the theory, namely 

Z = j VA^VA^VWVUVd'V^VdVe e-^^^'^'^^^^ (19) 

where 5 is given by ()18p . Now, taking the limit a — t- and integrating over the Lagrange multipliers 
{W, b^), one gets 

Z = J VA^VA^Vc^VcfVc'Vc' 5{D''^ Al)5{d^A'^) e-^YM+/d*x [c^m^^c^-c^ d,{d,c^+gr'^^Ay)] _ ^20) 
To deal with the diagonal ghosts (c*, c*) we perform the following change of variables 

^ f = + gr^'^AIc' , ^ f = , (21) 

with all other fields unchanged. Being linear in (c*, c?) this change of variables leads to a Jacobian which 
is field independent. Thus, we can verify that 

d^{d^e + gr^'^A^c^) ^ Cd^e ■ (22) 

Therefore, for the partition function we have 



j PfPf e-^'^*^^''^'^'^ j VA''VA^Vc''V(f5{DfA^^)5{d^A'^)e-^^^+^'^^''' 
N j VA^VA^Vc^m" 6{DfAl)d{df,A^^) g-^YM+J d'^x c-m-^c^ ^ (23) 
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with N being an irrelevant constant factor. Finally, integrating over the off-diagonal ghosts (c'^,c'^), it 
follows that 

Z = N j VA^VA' 5[DfA^^)5{d^,A'^) det(A^'^'') e"^™ . (24) 

As we shall see later, the expression ()24p will be taken as the starting point for the implementation of 
the restriction of the domain of integration to the Gribov region Vt. 



2 The Gribov region and some of its properties 

2.1 On the existence of Gribov copies in the maiximal Abelian gauge 

Before entering into the analysis of the Gribov region fi, let us spend a few words on the existence of 
copies in the maximal Abelian gauge. In order to deal with this issue, let us first take a field configuration 
(A^,^^) fulfilling the gauge conditions (jlOp . (jlip . Then, we can ask if there exists a gauge transformed 

configuration, (AJ^,^^), which fulfills the same gauge conditions as (^^, A^), i.e. 

D-^{A)Al = Q a^l;, = 0. (25) 

In the case of infinitesimal gauge transformations, this question can be answered in the affirmative if the 
following conditions can be verified for some value of the gauge parameters [lo"" ^u^): 

M''\A)u'> = , -d^{d^J + gr'^A^u') = , (26) 

where is the Faddeev- Popov operator given by (llSp . Notice that the diagonal parameter uj^{x) is 
completely determined in terms of the off-diagonal quantities uj"'{x) and 74^(x). In fact, from eqs. ([26]l . it 
follows 

= -ar^'^.d, [a^J) . (27) 

Thus, we can say that the condition for the existence of equivalent field configurations, or Gribov copies, 
relies on the existence of zero-modes of the Faddeev-Popov operator A4"'^{A). Explicit examples of such 
zero modes can be found in |31j . 



2.2 Determining the Gribov region Q, in the maximal Abelian gauge 

The so-called Gribov region Vt in the maximal Abelian gauge can be defined in a way similar to the 
case of the Landau gauge [3 [8] . More precisely, the region Q. consists of all field configurations which 
are relative minima of the auxiliary functional (jl2p . which is the same as to require that M"-'^ has to be 
positive, namely 

5'^F[A] = 2 J d'^xu^M''''{A)u^ > ^ M^^A) > 0. (28) 

Thus, we can write 

n = {^^, A^^ I D'l^'A'^ = 0, d^A^^ = 0, M'^'iA) > 0} . (29) 

The boundary 517 of Q is the region in field space where Ai"'^ achieve its first vanishing eigenvalue. This 
boundary is often called the first Gribov horizon, or simply the Gribov horizon. Notice also that for very 
small values of the coupling constant g, corresponding to the perturbative regime, the Faddeev-Popov 
operator behaves essentially like the 4-dimensional Laplacian, Ad""^ ~ —d'^S"'^, exhibiting in this case only 
positive eigenvalues. This means that perturbation theory is contained within the Gribov region O. 
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2.3 Some properties of the Gribov region 



In eq. (|29p . we defined the region Q as tlie set of fields fulfilling the maximal Abelian gauge conditions 
and for which the operator Al"^ is positive definite. Let us now establish some properties of this region. 
Let us first take a look at the Faddeev-Popov operator Ai""^, which can be written as 



M 



ab ^ab ^ab ^ab 



where 



Of = -D'^^Df , Of = gr^A'^^Df , Of = g^r^f^'A^Al . 



(30) 



(31) 



As pointed out in [3T1I17|, the operators Of and Of' are positive definite. We can easily show this by 
the following elementary calculation: 



\Oi\ 



> 0. 



(32) 



I as 



d^xir)^g^r'f'"''A';,Af,i^'= I d^x{gr'Air)\9f'''Al^'>) 



gr'^Air 



Concerning the operator Of, we note that its trace in color space vanishes, i.e. 



> 0. 

(33) 



tr02 = or = gr^A'^.D''; = gr''A';,i5'"'d^ - gf'^'A^^) = -gr'A'^d^ - g' r^A'^^A'^ = . (34) 

=0 

As a consequence, Of has both positive and negative eigenvalues. As we shall see in next sections, this 
feature can be used to establish a few properties of 



2.3.1 The region U is unbounded in all diagonal directions 

To prove this statement it is sufficient to observe that the purely diagonal field configuration (0, ^^) with 
A''^(x) transverse, d^A^^ = 0, fulfills the maximal Abelian gauge condition. Moreover, for this kind of 
configuration, the Faddeev-Popov operator Ai""^ reduces to the covariant Laplacian 

M''\0,A') = -D';fiA)D';^iA), (35) 

which is always positive for an arbitrary choice of the transverse diagonal configuration A^(x). We see 
thus that one can freely move along the diagonal directions in field space. The Faddeev-Popov operator 
will never become negative, meaning that 17 is unbounded in the diagonal directions. 



2.3.2 The region is bounded in all ofT-diagonal directions 

To establish this features, we observe that if {B^, BJ^) is a field configuration fulfilling the maximal Abelian 
conditions, D'fB^^ = 0, d^B^^ = 0, then the rescaled configuration {XB'j^, Bji), where A is an arbitrary 
positive constant A, obeys the same conditions, i.e. 

DfiXBl) = XDfBt = . (36) 
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Let now (A^,^^) be a configuration which belongs to $7, namely 

{^P\MiA'',A'M) = (VIOilV') - (V'lOslV') - ii^Psli^) > , (37) 

and let us evaluate the quantity {'^p\A4\^lJ) for the rescaled configuration (AA^,A^). From expressions 
(j3T]) one obtains 

{i;\M{XA'',A'm = (^IdlV') - Hi^\02\tP) - A'(V'|03|^) > . (38) 

Two cases have to be considered, (i) When < A < 1, the quantity above is always positive since 
{^l^\Ol\'^^J) > for every configuration and (A^,^^) belongs to the Gribov region by hypothesis, see 
eq. (f37|) . (ii) When A is larger than 1, we have to pay attention to the value of {^p\02\1p)■ If this quantity 
is positive, we just need to take A as large as we wish and certainly a zero-mode will be achieved. Notice in 
fact that the third term in eq. (|38p is always positive, ('i/'ICalV') > 0- If (^|C'2|V') turns out to be negative, 
it would be sufficient to observe that the contribution coming from {■ip\02\il^) is of first-order in A, while 
the contribution coming from (TplO^lip), which is always positive, is of the order A^. Thus, we conclude 
that even if ('(/'|02|^) is negative, we shall still achieve a zero-mode of 7W"* for A sufficiently large. In other 
words, moving along the off-diagonal directions parameterized by the rescaled configuration (Aj4^,Ap, 
with (A^, A^) belonging to the Gribov region Q, one always encounters a boundary dQ, i. e. the horizon, 
where the first vanishing eigenvalue of the Faddeev-Popov operator appears. Beyond 50, the operator 
Ai"''' ceases to be positive definite. 

Before ending this section, it is worth to discuss the case of a particular field configuratioiiH which 
can be handled by making use of the charge conjugation invariance displayed by the maximal Abelian 
gauge, see [331 SH ESI HTJ HSl |l9] for a detailed account on this symmetry. This field configuration 
would correspond to the case in which the eigenvalues of the operators O2 and O3 both vanish. If 
this particular configuration obeys the maximal Abelian gauge conditions p^ . pT]) . and allows for an 
arbitrary rescaling of the off-diagonal components of the gauge field, then it would give rise to off-diagonal 
unbounded directions in field space. Luckily, it turns out that this kind of configuration can be excluded 
by invoking the charge conjugation symmetry. In order to illustrate this issue, let us consider in detail 
the simpler example of SU{2). In this case, the operator O2 is absent and the eigenvalues of ^^{A) are 
given by the characteristic equation: 

g^AlAl - ss -g'A^Al 

-g^Aj^Al g'Al,Al-es 

For configurations of the form 

4 = ^4' (40) 
with arbitrary real A, it is easily seen that the characteristic equation displays in fact a zero eigenvalue. 
Furthermore, this configuration can be made to satisfy the maximal Abelian gauge conditions in a 
nontrivial way. It would appear then that by choosing A as large as we wish we could make the Gribov 
region Q unbounded along these particular off-diagonal directions. However, it turns out that Yang-Mills 
theory in the maximal Abelian gauge possesses a charge conjugation invariance, which is also enjoyed by 
the horizon term implementing the restriction to the region [431 [32| Wl\ HHj 139] . This symmetry is 
defined by 

4^4, 4^-4, 4^-4. (41) 

and is clearly not obeyed by the configuration (j40p . As such, these configurations can be excluded from 
the field space relevant for the Gribov region Q. Although the general case of SU{N) looks more involved, 

^We are grateful to Nicolas Wschebor for having pointed out to us this possibility. 



= (5^44 - e,){g'AlAl - .3) - 5^444^^ = • (39) 
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we expect that a similar reasoning applies as well. Needless to say, the charge conjugation can be also 
introduced for SU (N) . 

2.3.3 The convexity of the region O 

Let us now discuss the issue of the convexity of the region $7. Due to the nonlinearity of the gauge 
conditions, this property will be established for field configurations having the same diagonal components. 
Namely, consider two configurations (i?^,^^) and which obey the gauge conditions 

Df{A)B';, = 0, Df{A)C'^ = Q, 9^4 = 0, (42) 

and belong to the Gribov region, i.e. 

M^^B, A)>0, M^^C, A)>0. (43) 

In order to establish the convexity of the region $7, we should show that the field configuration {E'^, A^^) 
defined by 

i?^ = ai?^ + (l-a)C^, 0<a<l, (44) 

belongs to fi, i.e. M'^^{E,A) > 0. In fact, one can verify straightforwardly that 

7W"^(S, A) = a M^^B, .4) + (1 - a) M^^C, A) + a(l - a) Of{B - C) . (45) 

Since Of' is positive definite, as already mentioned, and since both (S^, and (C^, belong to 
we finally conclude that 

M''\E,A)>^, (46) 
showing that (-B^, belongs to $7 too. This proves the convexity of Vt. 

3 Gribov's statement about field configurations close to the horizon 

In [1], Gribov proved, in the case of Landau gauge, that for any field configuration located within the 
Gribov region and close to the Gribov horizon, there exists a nearby copy located on the other side of 
the horizon, i. e. outside of the Gribov region. This result was extended in [32] to the maximal Abelian 
gauge in SU{2). Here, we provide the necessary generalization to SU{N). First, let us consider a field 
configuration {C^,CJj^) lying on the horizon, that is 

<(C)C^ = 0, d^Cl = 0, (47) 

and 

M''\C)(l)l = , (48) 

where 4>q{x) stands for a normalizable zero mode of the Faddeev-Popov operator 7W"^(C). For later 
purposes, it turns out to be useful to introduce the diagonal components 0o(^) defined according to the 
second equation of (f26l) as 

^h = -gr''^d,ic;4). (49) 
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Let (A^, be a field configuration belonging to the Gribov region 0, and located close to the horizon 
dfl. For such a configuration, we can write 

^M = ^M + <' ^M = ^M + <' (50) 

where are treated as small perturbative components. As (^4^^,^^ obeys the gauge conditions 

(Iin]),(IIIl), from eq-dlZD it follows that 

<(C)< - gr'^a'^C'^ = , d^a^^ = 0. (51) 

The eigenvalue of the Faddeev-Popov operator corresponding to the field configuration (A^, A^^) is easily 
determined at first order in the small components (a^, a^) by using the standard perturbation theory of 
quantum mechanics, yielding 

+g^r'''f''a]^C'^4 - g^r'f'^Clal + a^C^)<A[5 ] . (52) 

The first term in the r.h.s. of (j52p vanishes because (C^, C^) lies on the Gribov horizon d^. 

We can now proceed as in [H [32], and introduce a third configuration (A^, A^), still located close to 
the Gribov horizon, defined by 

where 



= «;,-(9;.c/>U5r'*C'^'/'o), (54) 

are small as compared to (C^, C^). It is easy to verify that this new field configuration was constructed 
in such a way that it fulfills the maximal Abelian gauge conditions 

<(I)I^ = 0, 5^1;, = 0. (55) 

Therefore, we can ask ourselves if (A^,^^) is a Gribov copy of {A'^,A'^^), i.e. if these two field config- 
urations can be connected by a gauge transformation. In order to answer this question we shall repeat 
the argument of [U |32] , amounting to construct such a gauge transformation in an iterative way. Let us 
suppose thus that it is possible to connect the two configurations by a gauge transformation, so that 

Al = A^- {D-^J + gr^'A'^u:' + s/'^'M^') - \r'''u\Dfu:'' + gf^^'A^ + gf^Af^oj') 
-^r^^co^d^io' + gr'^'A'^^u^) + p-'^uj'iDl^'u;' + gf^'^A'^y + gf'^A'^y) + 0{u^) , 
^1 = id,u' + gr^'Ay) - ^r'^u-iD'^Lo' + gf'''A';^u'' + gf'^A<;^u^) + 0{lo^) , (56) 
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where, according to [T], the gauge transformation has been considered till the second ordeiH in the gauge 
parameters a;*,a;". Prom the maximal Abelian conditions it follows that 



-gr'^d^Lo^ + gr''^A'^y){D';uj' + gf'Ulujf + gf'^A 



0, 



0. 



(57) 



(58) 



The next step is to express {uj^,uj^) in terms of ((/>§, </)o). Let us start by writing 

uj'' = ro+ Vo . uj' = (t>l + vi, (59) 
with {ri°',rf) small when compared to Equation ()57p gives thus the following condition 

d^r^- = V''{C,a,<t^o) + Q-\C)rl>, (60) 

with T"^ and Q"* given by 



facdc r\ab 



2 racd fdbi i /~ic , 



2 fad fbdi / /~ic d , „C/^d\ib 



+Df{C) 



ybcd^c^jjde^cW^ ^ gfdefcy + gf-^c^M - ^f^cp'oidy + af^'cy,) 



9 fbci. 



+gf-icy) - g^^idy + gycy){D';{C)r, + gf'^C'^4 + gf^'C^) , (61) 



Qab ^ ^gfab^ci^Q^_ gjabc^cj^d, _ g2fac^fcbJci^C|^_ 

We notice now that equation (|6Up can be solved for 77" (x) in an iterative way: 



(62) 



(63) 



allowing us to obtain a recursive expression for the parameters (ct;",a;*) as well as for the gauge trans- 
formations (j56p in terms of ((/>g,(/>g). We have shown thus that the configuration (A^,^^) is indeed a 
Gribov copy of (^^, Ap. Furthermore, equation ([60]) can be used to establish another relevant property. 
Recalling that M"'''{C)(I)q = and that M""^ is Hermitian, we have 



(64) 



^Notice that if we had taken only the first order gauge transformations, we would get 



AO- *a / T\ab b , j-abc Ab c , rabi *b 



Al - {d^u' + gt'^AlLo') . 



Then, applying the maximal Abelian gauge conditions we would obtain the following equations: 

M^\A)u,' = , d^{d^u>' + gr'-A^co') = , 
which have no solutions since, by hypothesis, (j4|i,ylp) is not located on the horizon dil. 



12 



Therefore, from ()60p . it follows 

j4 la J ri „ fabi i T\bc//~(\ic „f(icd„c r~\ab / ib , 2 racd fdbi i r^c ib 2 fad fbdi / /~ic d , „c r~id\ 



~<dij^ 



^y_fbci^j^cd^C^^d ^ gfcde^d^e ^ gfcdjcf^^^^ 



(65) 



Moreover, the eigenvalue of the Faddeev-Popov operator corresponding to the configuration (A^, can 
be obtained in the same way as done before for (A^, j4^), amounting in fact to replace by a^. Thus, 



e{A) 



acd~c r\db/ 



+g^rT'alCl<Pl - g^r'niC^al + a^C^^;^, 



~id\ lb 



Finally, from equations 



and dMD we obtain that 

e{A) = -e{A)^ 



(66) 



(67) 



Since e{A) > by hypothesis, it follows that e{A) = —e{A) < 0, showing that the Gribov copy {A^,A^^) 
is in fact located outside of fl. 



4 Restriction of the domain of integration to the Gribov region Q 

As we have seen before, field configurations located inside and close to the boundary dO, have copies 
outside of 0. Therefore, restricting the domain of integration in the Feynman integral to the Gribov 
region allows us to get rid of a certain amount of copiec|. To implement this restriction we shall follow 
[H |32] and introduce in the partition function of the theory, eq. (|24p . the factor V(0) which formally 
constrains the domain of integration to the region 

2- = j dfie~^^^ = j dfie-^^^Vin) , (68) 

where dfi is the Faddeev-Popov functional measure given in eq. (j24p 

dfi = VA°''VA'''^^ 6{D ■ A°'')6{d ■ A'''"^)(det M) . (69) 

In order to obtain an explicit expression for V(il), we make use of so called no-pole condition [l], which 
is a condition on the connected two-point function of the off-diagonal Faddeev-Popov ghost fields. As 
pointed out in [l], see also [50] for a pedagogical introduction, the no-pole condition stems from the 
positivity of the Faddeev-Popov operator, Ai"''' > 0, within the Gribov region il, according to eq. (l29p . 

^Nowadays, it is known that the restriction to Q does not enable us to get rid of all copies. As already remarked, a 
further restriction to a smaller region, known as the fundamental modular region, should be implemented 12]. Though, till 
now, this seems to be beyond our present capabilities. 
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As a consequence, within Q, the operator A4"'^ is invertible. Moreover, its inverse, {Ai"'^)^^, which is 
nothing but the connected two-point off-diagonal ghost function evaluated in an external gauge field 
background: 

^ ''^ J PgofT-p^off dizcb{z)Mi"'{z)c'^{z) ' ^'^> 

is positive too within Q., i.e. {M""^)^^ > 0. Following [I], the next step to achieve the factor V(J1) is that 
of considering the quantity 

g[k- A) = j <fx dS G{x, y- A) e'^^'^-y^ , (71) 
where Q{x,y]A) stands for the trace, in color space, of namely 

> - N{N - 1) ~ iV(iV- 1) fVc'^^Vc"^ Jd^zc>'{z)M^Hz)c-{z) • 
Parametrizing thus Q{k]A) as in p!| 132] 

where J is constant, the no-pole condition can be stated as [1]: 

a{k-A)<l. (74) 

From this condition, it follows that the off-diagonal ghost propagator has no poles for any finite value of 
the momentum k, so it cannot change sign by varying k. It will always remain positive, meaning that 
the Gribov region Q will be never crossed. The only allowed pole is at A; = 0, where expression (j73p 
becomes singular, meaning that we are approaching the boundary where is in fact singular, 

due to the zero modes of M""^ . Moreover, taking into account that a{k; A) is a decreasing function of the 
momentum k [T], for the final form of the no-pole condition we might take [1] 

a(0;A)<l, (75) 

which is very suitable for explicit calculations. Therefore, the factor V{Q) implementing the restriction 
of the domain of integration in the functional integral to the region Vt is 

V{^) = e{l-a{Q-A)) , (76) 

where 6 stands for the step function. Consequently, for the partition function we can write 

Z = f d/ie-^YM= [ dfie-^^'^eil-aiO^A)). (77) 
4.1 Evaluation of V(fi) 

Let us face now the characterization of the factor V(n). To that purpose, we start with expression (j72p 
and evaluate g{x, y; A) order-by-order in perturbation theory. Performing the calculation until the third 
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or 



dein we obtain!!! 



Goix -y)-g^J <fz [N A'^{z)A\Xz) - Al{z)Al{z)] Go{z - y)go{x - z) 



xd'^'Goix - zi) go{z2 - y) dt'Qoizi - z^) 



+ 



N{N -1) 



d''z^d''z2r'''f'"''r'''AUz,)Al{z2)Al{z2) 



{d'tlG(){x - zi)^Qq{z2 - y)Go{zi - Z2) + Goix - Z2) Qq{zi - y) 9^'^o(z2 - ^1) 
d''zrd''z2d^zJsr''f'^r^AUz^)Ai{z2)Al{z^) 



where 



N{N -I) 

'd^^Goix - zi))^o(^2 - y)(dl^Go{zi - z^))dl^Go{z^ - ^2) 

d^q 1 



Gq{x - y) 



-iq{x-y) 



(27r)4 g2 



(78) 



(79) 



In momentum space we can write 



Gik;A) 



d'^xd'^y G{x, y; A)e 



ik{x-y) 



1 

+9' 
+ 



fc4 



d^q fAl^{q)Al^{-q) _ A-{q)A-{-q) 

N{N -1) 



N - 1 



d'^q kf,{k - q)y 
(27r)4 k^{k-qf 

f d'^p d'^q 



N{N -1) 

,3 



(iV-1) ^ 
ik 



Al{q)A}^{-q) + 



N 



-Al{q)Al{-q) 



+ 



N{N-l) ^ 

^abc rbde fadi ac 



+ 
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(27r)4(27r)4 VA;^(p+A:)2 k'^^k - p)^ 
d'^p d'^q ikfj{k — p) ^{k — p — q) 



r''nnAi{-p)Ai{-q)Ai{p + q) 



8r''f''^r''A^J-p)Aii-q)A';ip + 



N{N -I) J (27r)4 {2-kY k^{k -pY{k-p- qY 
+^r'^f'''f'''A',{-p)AiA-q)K{P + 9) + ^r'''f^'f'''Al{-p)Al{-q)A^,{p + < 
+4/'^''7'"7''^^'^^(-p)4(-'Z)4(P + 9) + 2r^^f^'f'''A'^{-p)Al{-q)Al{p + . 

+2/'^''7'"7'"^^^(-p)4(-9)^^(p + q) + 2r''f'''f'''''A';,{-p)Ai{-q)Ai{p + , 



(80) 



In [TJ[3^ the calculations were performed only till the second order, as this is sufficient to determine how the restriction 
to affects the form of the tree level propagators. Furthermore, a third order calculation will be needed in order to 
characterize the horizon function, as recently done in [51] in the case of the Landau gauge. 
^More details of this lenghty calculation can be found in the Appendix [51 
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According to eq. (l73p . the expression above can be rewritten as 

g{k;A) = ^{l+a{k;A)) + ^, (81) 
where the form factor a{k;A) and J7 are given by 

f d'^p d-^q ( ik^ - P) \ ^abc ^bde ^adi 



+^r''f'"''f'"''A^,{-p)AU-q)Al{p + q) + Ar^'^f'''f'^^A'^{-p)Al{-q)A^,{p + q) 
+4/'^'7'"7''^^'^^(-p)^U-'?)4(P + q) + 2r^'f'''f'''A'^{-p)Al{-q)Al{p + q) 
+2r''f'''f--A'^{-p)A'^{-q)Al{p + q) + 2r^-f''-f'^^Al{-p)Al{-q)^^{p + q) 

+r'"f'''f''fAl{-p)Al{-q)Ai{p + q)) , (82) 

, f dS fA^{q)A]A-q) A-{q)A-{-q) \ 
^-~^Jj2^^[ N-1 iV(iV-l) )■ ^^^> 



Notice that is independent from the external momentum k. Also, expression (j82p can be simplified by 
making use of the maximal Abelian gauge conditions in momentum space (jl3p . namely 



q,A^^ = , q,A-{q) = tgP'^ J ^ A^^{q - p)A'^{p) . (84) 



Then, for a{k; A) we obtain 



Af(JV-l) fcz; (2i)*(2ir)''(t-(j)2' "' 



+4/"'7'"V'"^4(-p)^i(-'Z)^^(p + q) + Ar''f'''f'''Al{-p)Al{-q)Ai^{p + q) 



+4/"'7'"7''^^'^^(-p)^U-'Z)4(P + <!) + 2r''f'''f--A^^{-p)Al{-q)At{p + , 

+2/"'7'"7''^^^^(-p)^U-9)^^(p + g) + 2/'^^7'"7"'^^^^(-p)^^(-g)^l(p + ^ 



+/"'7''7"'^^^^(-p)^^(-'?)4(P + 9)) + 0(/) . (85) 
Using the relation 

N — 9 

jahc jbde jadi ^ jcei (86) 



16 



and performing the change of variables p 
(j85]l . we get 



a{k;A)=g 
where 



Iff' 



.(1) . ig^{N -2) ^ , 



q and g — > p in the third hne of the r.h.s. of expression 



A;2 



(87) 



Jt^Jt^^ \{k - q)^ {k + qf 
d'^p d'^q 



N{N - r 

r^A';^{q)Ai{p-q)Al{-p), 

(2.)^ (2.)^ (. -pnk-p- q? y^r''f''^r'^i-p)^(-^^(^ + ^) 

+4/'^''7'"^/"'^^^;(-p)^i(-g)^^(p + + ^r''f'"''f''''A^,{-p)Al{-q)Aiip + < 
+4/'^''=/'"^7"'^^'^^(-p)4(-9)4(P + 9) + 2r'^f'''f''-A'^{-p)Al{-q)Al{p + , 

^2/afec^6d^^dae^c (_^)^i^(_^)^e + + 2/'^'''^ f'^' f'^' A';^{-p)Al{-q)Al{p + 

+4/'^''7'"7""^'^^(-p)^U-9)4(P + 9) + 2r'^f'''^f''''-A'^i-p)Ali-q)AUp + q) 
+2r''f'''f'''Al{-p)A',{-q)Al{p + q)+ 2r''f'''f'^'Al{-p)Al{-q)A',{p + q) 



(89) 



(90) 



+r'''f'''f''fAl{-p)Al{-q)Ai{p + , 

In order to be able to evaluate the form factor a{k;A) at zero momentum, eq. (|75p . 

a{^;A) = lima(A;;^). 



(91) 



(92) 



we need to analyze the limit k ^ of the integrals (I88p - (I9ip . As these integrals depend explicitly on the 
Fourier components of the gauge fields, 74^(g), j4^((jf), the existence of the zero momentum limit. A; — ?> 0, 
of expressions (|88p - ()9ip relies on the infrared behavior of the fields A^^(q) and A^(g). For example, 
considering (j88p . it is apparent that the zero momentum integral 



d^q 1 



{2tt)^ q^ \{N - I) 



does exist provided the quantities A^^{q)A'^^ 

Al{q)Al{-q) 



-^U9)A 



N 



-M{q)Al{ 



N{N - 1)"^' 
A^(g)^^(— (7), are not singular a.i q'^Q, i.e. 



(93) 



2\a 



a> -2 , 
/3>-2, 



(94) 



so that expression (|93p is integrable at g ~ 0. At the present stage, a formal proof of the infrared behavior 
expressed by (j94p cannot be given. An explicit check of the validity of ()94p can be done only after the 
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evaluation of the gluon propagators. In other words, as done in the original paper by Gribov [Ij, ()94p is 
assumed as working hypothesis in order to ensure the existence of the zero momentum limit of expressions 
([88]) - (fOTI) . Afterwards, this assumption has to be checked out by looking at the behavior of the diagonal 
and off-diagonal gluon propagators at the origin. Nevertheless, we can provide a physical justification 
for assuming (194p . which follows from the implementation of the restriction to the Gribov region ^1. As 
also observed in [l] in the case of the Landau gauge, the restriction to the region Q, amounts to put a 
boundary in field space, namely the Gribov horizon, which manifests itself through the appearance of a 
dynamical mass parameter 7, which is known as the Gribov mass As we shall see in the next sections, 
such a parameter will provide a natural infrared scale for the field configurations, resulting in a gluon 
propagator which is not singular at the origin, thus justifying the starting working hypothesis ()94p . In 
fact, looking at the expressions for the diagonal and off-diagonal propagators given in eqs. (I209p . (|210p . 
one sees that they are not singular at the origin. 

The integral lj^\k) is automatically vanishing when k = 0; the integrals /^l/''(0) and must be 

proportional to 6^1, in the limit /c — t- 0, since 5^^, is the only Lorentz invariant second rank tensor avail- 
able; and finally iju/aik) must vanish when A; — >■ 0, as in four dimensions there is no Lorentz invariant 
third order rank tensor. Summarizing, we have 



(0) = ^ 5^u , (by Lorentz covariance) , 

(2) 

III {0) = , (automatically) , 

I/iua (0) = , (by Lorentz covariance) , 

ij^) (0) = ^'^^ , (by Lorentz covariance) . 



(95) 



Thus, taking the limit A; — )■ of a{k; A), for the form factor a{0; A) we get 
a{0;A) = limcT(A:;^) 

fc-i-O 

= ,.il>)__«L_Z>)+o(/) 

^ 4 N{N - 1) 4 ^ ^ 
2 f d^Q ( 1 A\{q)A^^{-q) N-2 A'i{q)Al{-q) 



(27r)4 V(A^ - 1) g2 4N{N-1) q^ 

ig^ f d'^p d'^q 



NiN - 1) J i2.f i2.f p^ip + q? yy'''f'^r^^^i-v)Ali-,)A\^ + .) 
+r''f"''f'''M-p)K{-q)Al{p + q) + r^-^f'^-fd-^A'^{-p)A'^{-q)Ai{p + q) 
+r''f''f'''^i{-p)^',{-q)A{{p + q) + \r'''f'"'f'''A^x{-p)A'^{-q)Al{p + q) 
+ lr'"f'"''f'^'Al{-p)A'^{-q)Al{p + q) + ^f->'-f<i-f'^^Al{-p)Al{-q)A\{p + q 



+ lr''f''''f''^Al{-p)Al{-q)A{{p + q)^ . (96) 
As a useful check, we notice here that, for N = 2, the previous results obtained in [32j are recovered. 
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4.2 Gribov's quadratic approximation and the gap equation 



Having characterized the factor a{0;A), eq. (|96p . we start to analyze the consequences on the theory of 
the restriction of the domain of integration to the Gribov region, as stated by eq. (l77p . As we shall see, 
this will amount to modify the starting gauge field theory by the addition of a nonperturbative term, 
usually referred as to the horizon function [71 [U |48]. To illustrate this point, we discuss first the form 
factor a{0;A) in the so called quadratic approximantion i.e. we shall first consider a(0;A) till the 
second order in the gauge fields. Thus 



a{0;A) 



^ V^\(N-l) 

q \y I 



1 A\{q)A\{-q) , TV -2 Al{q)Al{-q) 



+ 



4Af(Af - 1) 



(97) 



where, according to [T], a finite volume V has been considered. As in the case of SU (2) [32], the quadratic 
approximation will enable us to obtain the first term of the horizon function, which will play an important 
role for its all order extension. 



Using the integral representation of the step function. 



for the partition function (j77|) we get the following expression 



27riC 



d[i e 



-S'ym+C{1-'^(0;A)) 



(98) 



(99) 



Moreover, in the quadratic approximation, Z — > Z^^^^: 



-'quad 



^^.(fc;C,«) 
Q;f.(^;C,/3) 



lim 

«->o / 21X1(1 



+ 2iV(iv-i)y/fc2 'V 



1 I knkij 
a ' 



1 ^ ) k^ki, 



{N-l)Vk^J 
where use has been made of the relations: 

6{D-A°^) oc lime-2^/'^'^(^"^M)% 

a— >0 

did-A"'-^) oc lime-i^'^*"^^'^^^)', 

13^0 

^5]A;l■^(A:)e-'=^ 



(100) 



A'l^\x) 

^k,k' 



y Jv 



(101) 
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Integrating over the gauge fields gives 



-'quad 



lim 



dC 



/(C;«,/3) 



C-lnC 



1 

2 

3N{N-1) 



^ _ In C - ^ In det 7^2^ (A:; C, a) - - In det Q^, (A:; C, /3) 



E 



In A;^ + 



3{N - 1) 



E 



ln( A;2 + 



(iV- l)yA:2 



2iV(iV - l)yA:2y 3 \a 2N{N -l)Vk'^ J 

(102) 



1 

+ - In — + 



/3 {N-l)Vk'^ 



The integral over can be evaluated in a saddle point approximation [T] 

Z,.., ^ lim e/(C*^"'/^) , ^-^^^^ 



a->0 



0. 



Condition (|103p gives us the following relatior0 

39' 



1----E 



1 ^ 3(iV - 2)^2/4 



-y 



y , 4 , 2.g2C* y , 4 {N-2)gH* 

^ fc (Ar-l)y fc 2N(N-1)V 







(103) 



(104) 



(Ar-l)y k ^ 2N(N-1)V 

As in [T], in order to take the thermodynamic limit, V — )• 00, we introduce the so called Gribov parameter 
7 

(105) 



4 r C 
7 = hm 



Therefore, in the infinite volume limit, V ^ 00, we get 



+ 



d^k 3(7V- 2)^2/4 



(2-)^ k^ + '4^' J (2vr)4 fc4 + (^.^ 



(106) 



where we have neglected the term 1/C*- As observed in [T], the Gribov parameter is not an independent 
parameter of the theory. It turns out to be determined at the quantum level in a self-consistent way 
as a function of the coupling constant. In fact, eq. (1106p is nothing but the gap equation enabling us to 
express 7 as a function of the coupling constant g and of the invariant scale Aqqd, see ref.[l3]. 

Coming back to the partition function in the saddle point approximation, taking the form factor a{0; k) 
in the quadratic approximation, and performing the thermodynamic limit, we can write 



Z = M' j dfie" 



(Sym+Sg) 



(107) 



Here, TV' is an irrelevant normalization factor and is a nonlocal term, given by 



Sg = lim Ca{0;A) 
= / d'xd'y 



2 . TV — 2 

AUx)go{x - y)AUy) + Al{x)go{x - y)A-(y) 



N-1 



/ d-x 



4* 

N-1 



1 



-92 



N -2 
^ 2N{N - 1) 



2N{N - 1) 
1 



-52 



(108) 



'^Notice that the gauge parameters a,/3 have been set to zero, see ea. p01[l . 
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Equation ()107p shows that the implementation of the restriction of the domain of integration to the 
Gribov region fl can be achieved by adding to the Yang-MiUs action a nonlocal term, given in eq. ()108p . 
Expression (jlOSp represents only the first term of what is usually called the horizon function, whose final 
expression will be discussed in details in the next section. 



5 Characterization of the horizon function 



5.1 A few words on the horizon function in the landau gauge and in SU{2) maximal 
Abelian gauge 

Before entering into the technical details of the determination of the horizon function for SU{N), it is 
useful to give a short survey of what has been already achieved in other cases. In the Landau gauge, it 
has been shown by Zwanziger [71 [8] that the horizon function implementing the restriction to the Gribov 
region jil is 

^Landau =9^1 d^X A Df{x) [{-d ■ (x, (y) , (110) 

where [{—d ■ D)~^]^^ {x,y) is the inverse of the Faddeev-Popov operator 

- id . D)^^ = -d,Df = -d,{d,5^^ - . (Ill) 

Accordingly, for the partition function implementing the restriction of the domain of integration to the 
Gribov region, we write 

^Landau = / e"^™ = / d/, e-(^™+^'^"") , (112) 

Jn J 

where, to the first order, the Gribov parameter 7 is defined by the gap equation [TJ [H] 

^ ~ y (2^ FT2]vy • 

One should observe that expression (jllOp is nonlocal. Willing to give a quantum field theory meaning to 
such an expression, one should be able to cast it into a local and renormalizable form. We mention here 
that a consistent, local and renormalizable framework for expression (jllOp is in fact available, see [11] for 
an updated discussion of this nontrivial issue. Moreover, it is worth to point out that Zwanziger horizon 
function (jllOp turns out to be equivalent to Gribov's ghost form factor, as shown in [51] by comparing 
the expansion of both expressions till the third in the gauge field. 



In the case of the SU{2) Yang-Mills theory in the maximal Abelian gauge, the results found in 
are very similar to the ones obtained by Zwanziger in the Landau gauge. Summarizing, the form of the 
horizon function obtained by requiring that it can be cast in local form, that it can be renormalizable 
and that it reduces to Gribov's ghost form factor in the quadratic approximation, is given by 

Hmag=9'I d^xd''xe'''A,{x){M-'r'{x,y)e''A,{y), (114) 



*We remind here that, analogously to the case of the maximal Abelian gauge, the Gribov region in the Landau gauge is 
defined through the positivity of the Faddeev-Popov operator, i. e. 

n^{A^, - d,A^ = 0; - gf^^^Af,) > } . (f09) 
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where {M ^)°'^{x,y) is the inverse of the Faddeev-Popov operator 

^ab ^ -DifDf - g^e^'^e^'^A'l^A'l . (115) 

5.2 Candidates for the horizon function in SU{N) 

Let us face now the problem of finding the horizon term for the general case of SU {N) . The first step 
for accomplishing this task is that of selecting the possible candidates which are compatible with the 
following requirements: 

1. it should give back to the Gribov term (jlOSp in the quadratic approximation; 

2. it should reproduce the particular case ()114|) when N = 2; 

3. it should be localizable and compatible with the requirement of power counting renormalizability. 

Five possible candidates which are compatible with the aforementioned requirements have been found, 
namely 

= N^frT) I f'^l,^i^~')tyf''''K,y + 2N{N-l) / (H^) 





1) 






N{N - 


1) 






N{N - 


1) 






N{N - 


1) 



+{N-2)Cg^f r'=M^,.(X-^)SV^<,, (119) 

J X,V 



rSth 



2N{N - 1) 



x,y 



In all expressions (fTTHI) - (fT20]) . the nonlocal operator {M~^)f^y is the inverse of the Faddeev-Popov 
operator (jlSp and the symbol ^ stands for f d'^xd'^y. In eq. (jll9p . the parameters ^ and C are positive 
and obey the following constraint: 

In expression ()120p . the parameter /3 is positive and is related to a according to 

a'^ + {N -2)(3 = N -2. (122) 
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Let us now see how these five candidates fulfill the three requirements. The first criterion claims that, in 
the quadratic approximation, the horizon function has to reduce to the Gribov term defined in (llOSp . In 
the quadratic approximation, the inverse of the Faddeev-Popov can be taken as 

^M-'r'^^. (123) 
Moreover, using the following relations between the structure constants 

jabi jabc q j^abi jabj jV(5*J jabc jabd ^jy 2)(^'^'^ jaic jaid ^cd (124) 

one can easily check that all five candidates (jll6p - (|120p obey the first requirement. 

The second requirement states that, when N = 2, the horizon function must equal expression (jll4p . 
It is easy to check that this requirement is also fulfilled by observing that in the case of SU{2) we have: 

jaM _^ jab3 ^ ^ab ^ jabc ^ q ^ ^ = . (125) 

Besides, there are terms proportional to [N — 2) which automatically vanish for N = 2. 

The third requirement is about the possibility of the nonlocal horizon terms be cast in a local form 
by introducing a suitable set of auxiliary fields, in a way compatible with power counting renormalizabil- 
ity. Let us analyze each case separately. The first candidate (jll6p consists of the sum of two different 
terms which can be localized independently as follows: 

X j VXDX (det M)4A^(^-i) exp - j (fx {xi'W^X^^ + y /"^"A;^(A - A);^''^ 
where O"^ stands for 



, (126) 



(127) 



N{N - 1) 



In this expression and (A^^,AJ^^) are two pairs of complex commuting auxiliary fields. The 



determinant, (det A^)4^(^-i) , can be exponentiated by employing a pair of complex anticommuting 
fields (c^^^O: 

(detX)4^(^-i) = J P^PcDe/'^'--S^-^"''< . (128) 
Thus, the local version of H^^^ is 



-f^Local — I d'^X 



.T^f^c \ /fab, he -ac \/tah, he , n2 rabi At /,„ ,-\c 



+Xl'M'^'X'^' - r^'M-'ii"; + ^ r'^A'^^iX - A);i^] , (129) 

where , fj'j^) is another pair of complex anticommuting fields similar to ((^^^, w^^). As one can easily 
see, the set of fields {ip,(p,u},Ld) is responsible for the localization of the diagonal sector, while the set 
(A, A, T], f)) allows us to localize the off-diagonal sector. 
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In a similar way, we can localize the remaining candidates. To localize expressions ()117p and (jllSh 
it is necessary to employ only one set of localizing fields, (99, w, a)): 



-Hi 



2nd 
Local 



r_r3rd 
-"Local 



(fx 



—ac KAab he ,-,ac \^ab, he , q2 I rab% ai , ^ fabc Ae \ ,^\a.l 



1 



?abc Ae 



\ab 



ah 



(130) 
(131) 



In the case of expression (jll9p we find 



4th 
Local 



,^fic k^ab,bc -ae KAab, be , q2 rabi Ai / , ,^\ab 



ab 



(132) 



Notice that, here, a new set of fields has been employed: {ip'^)''P'^) stands for a pair of complex com- 
muting fields and {x^^X'^^) for ^ P^ii^ of complex anticommuting ones. We also observe that these fields 
carry both off-diagonal and diagonal indices. 

Finally, the local version of (|120p is 

(fx 



rrSth 
-"Local 



Cjoej^ab^be ^ , ^ab^ j^. ^ ^r^A^ ] - ^) 



cabi At 



a 



abc Ae 



\ab 



bi 



(133) 



It is worth to point out that all local expressions, eqs. (|129p . (jl30p . (|13ip . ()132p . ()133p . are power counting 
renormalizable. All five candidates, (|116p - ()120p . fulfill the three requirements given at the beginning of 
this section. 



A further requirement is thus needed in order to select only one candidate. The natural set up is that of 
strengthening the requirement of the equivalence between Gribov's form factor a{0; A) and the possible 
horizon terms, as done in [51J in the case of the Landau gauge. This amounts to expand the horizon 
terms, eqs. (|116p - p2Up . till the third order in the gauge field and compare the resulting expressions with 
that which we have already obtained for a{0;A), see equation ((96]) . 



5.3 Selecting only one candidate for the horizon function 

In this section we show that the only candidate which is compatible with the expression of cj(0; A) is 
^2nd^ eq. (|117p . To that purpose, the expansion till the first order in the gauge field of the inverse, 
{A4~^)°'^, of the Faddeev-Popov operator 

>1"''(x)(M-i)^'=(x, y) = 5^-5{x - y) , (134) 

is needed. Recalling that M^-^^x) is defined by expression ([15]), we can write 

^ab ^ _^abQ2 ^ 2gr^A^^df, + gr^'A^^df, + 0{g^) . (135) 
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To evaluate the inverse of Ai""^ order-by-order, we set 

iM-'r\x, y) = gS\x -y)+g Qf{x, y) + 0{g^) , (136) 
so that, to the first order, we have to solve the equation 

M^'ix) {g^\x -y)+g gf{x, y) + 0{g^)) = 6-'6{x - y) , (137) 
where Qq'{x — y) is the free ghost propagator: 

5SV - !/) = = / pifi ? ■ (138) 

To the first order, we have 

- dig fix, y) = -2nAl{x) dig^\x -y)- r''Af^{x) d;g^\x - y) , (139) 

which gives 

QO,b ( ^ „,^ o / j4 ^ ^ I fabi A'i ( ^\ , ^ fabc /ic / " ' ^ 

Therefore, 



gr{x,y) = -2 / rf'^^3^^r"4W + 2^" ^^M(^)J^M^3^- (140) 



iM-^T\^^ y) = -2j d'z j—-^ i^r'^A^^iz) + -r^^^^(^) j ^M^f^ + 0(9') ■ (141) 

We are now ready to expand the horizon function (1117P in powers of the gauge field. Introducing the 
combination 

one can write 

H'^^ = ^) / d'xdSBl^{x){M-^rHx,y)Bl\y) , (143) 

which, by means of eqs. p4ip and (|142p . becomes 

[ d^xd'^yd^z Bl\x) Bf{z) dly^^ B^^iy) + 0{g^) . (144) 

J \x z\ \z y\ 
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N{N - 1) 
Explicitly, we have: 

1 I f I- - 2 o^T ^T — ^\ / ^ U. _ i, 2 



J d^xdSd^zj-^d^^j^^ (/-7"'v'"='4(^)^i(^)4(y) 



N-lJ 2N{N-1)J |x - 

4g^ /■ j4„.j4..^4., 1 1 

iV(A^ - 1) 

+lnr'^f"'^A'^{x)AUz)Af,iy) + ir-r''7'^H(^)^^(^)4(y) 

+\r'r''f''^A^^{x)Al{z)Al{y) + i/-'^/'^^V''^'<(x)^i(^)4(y) 

^l_^acd^a6i^bc/^d(^)^, + ^rT"' f'"' Af^{x)Al{z)A>;{y) 

+lr''r''V''^Af,ix)Aliz)Aliy)^ + 0(/) . (145) 
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Taking the Fourier transform 



2g^ f d\ A\{q)A\{-q) , {N-2)g'' f d\ Al{q)Al{-q) 



(N-l) J (27r)4 g2 2iV(iV -1)7 (27r)4 g2 

+r''f"''f'''''A^,i-p)Ai{-q)A';,{p + q) + r^^f'^-f-lA^^{-p)Al{-q)A^^{p + q) 
+r''f''''f''''^A';,i-p)Ali-q)Aj,ip + q) + }-f-'^f'^-f-^Al{-p)Al{-q)Al{p + q) 

+ \r''f'''f'''Al{-p)A',{-q)Al{p + q) + lrb^f'^-f--A'^{-p)Al{-q)Al{p + q) 

+ \r''f'''f''^Al{-p)Al{-q)Al{p + g)) + 0(/) , (146) 

and comparing with Gribov's form factor o"(0; A), eq. (l96p . it is apparent that 



cj(0,yl) = ^^H"^""^ + 0{A^). 



(147) 



The same procedure can be easily repeated for the other expressions, H^^^ , H^^^ , H^^^ , H^^^ . The cor- 
responding expansions do not coincide with equation (j96p . Therefore, from now on, we shall consider 
expression (jl43p as the correct horizon function we were looking for. Consequently, for the partition 
function defining the theory we shall take 

Z = [ d/.e-(^-M+7^//-'^). (148) 



In the next section, we shall construct a complete local action out of i^^nd shall study the possible 

dimension two operators which can be introduced. The effect of the condensation of these operators on 
the tree level gluon and ghost propagators will be also worked out. 



6 Local dimension two operators and the tree level gluon and ghost 
propagators 

6.1 A local and invariant action 

According to expression ()148p . the general local action we can construct for the SU{N) Euclidean Yang- 
Mills theory in the maximal Abelian gauge reads 

Z = j e~^° , 

Sq = 5yM + 5'mAG + -SLocal j (149) 

where 

p$ = VA°'VA'Vb°'Vb'Vc''V^Vc'Vc^VipVipVuVui , (150) 
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and 5Locai 
So = 



-^Locai given by eq. (|13Up . Explicitly, we have: 



eabi a a h\ 



-ac A ^ab, he 



+ 



27^ 



N{N - 1) 



eabi ai ,9 rabc 



ab -ab 



(151) 



Having at our disposal a local action, we may investigate its symmetry content. Let us start with the 
BRST symmetry, already obtained for the gauge fixed theory, see eqs. (fT6]) . As established in the case 
of the Landau gauge [71 13 [13] and SU(2) maximal Abelian gauge [l9], the auxiliary fields {ip,ip,uj,Cj) 
transform BRST quartet, i.e. 





— 1 




= 




— ^fl ) 




^< 


= ^fi > 




= 




= . 







(152) 



As a consequence, we get 



he 



-ae A rah , he 



-ac \ fab, ,be\ 



s I d^x Bfi^f - ^l') = 9' I d^x [{sB^;){ipl' - O + i?^''^] 



^ab\ 



ab\/ab 



jab, nbi 



(153) 



where 9"^ is given in eq. ()127p . B'^^{x) in eq. ()142p . and the operator A/''^'' is defined as 



J\^ab^b A Aab^b\ \ ^ab 



2gr^\sA^^)d^ + gr''{sA'^^)d, + gr'^d.isA^^) + 2g^ f'\sA^^)A^^ 



,2 fade fhdi „/■ Ae At \ „2(facifbdi , radi fbci\ 



- '2gr\dy + gr'''A'^y)Df - gr^D^'c' + grU^cf + gnAld)D\ 
+gr''d,{dy + gr^'A'^) + 5V"'7''M^(a^c^ + gff'Alcf) 



db 



_g2^jaeijhdi ^ pdijhei^^jjee^e ^ gf^efj^e^^f + gf^^JAl(^)A 

with representing all the off-diagonal fields of the theory 

r /la ua ji -a ,ab ,-ab , ah -afei 



(154) 



(155) 



The first term S^^^ can be rewritten in a BRST invariant way by making a linear shift in the variable 

<(x) [a m [111 149], 

<''(x) ^<^(x)+ / A(X-')"'^(x,y)AA^''(y)95f(y). (156) 
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Thus, 

d^xCj^^'^M^^ipl^. (157) 

In order to deal with the second term S^'^\ let us notice that 

B^'^f = Bfsof = s{Bfuf) - [sBf) < , (158) 

then, 

5(2) =e^ j d'x [Bf^f - siBfu^f ) + {sBf) u^f] . (159) 
The last term of the r.h.s. of (|159p can be eliminated by performing a second linear shift in uj'^{x) 

[3 El m 119] 

ojfix) ^ ojfix) + e'J d'y {M-'r{x, y) sBf{y) . (160) 
Therefore, the local version of the horizon function might be taken as 

^Locai = s j d^xO^^'^M'^V,' + 0'J d'x [Bf^f - .(«)] . (161) 

Consequently, the action Sq becomes 

5o = 1 d'^x ( + \f;,FI^ + ih'^ OfA^ - ^M'^'c' - {DfAD [gr^'^c^c'^ + gr^'c^d 

+ sjd^x O'lfM^V,' + e'J d'x [Bf^f - s(«)] . (162) 

As it happens in the case of the Landau gauge [71[8l[13] and SU (2) maximal Abelian gauge [l9], the action 
does not possess exact BRST invariance, which turns out to be broken by soft terms proportional to 
the Gribov parameter, i.e. to 9"^. In fact 

sSq = ff' l\\,reak i (163) 

where 

« - {^""""'{d^c'+gf^Alc^) + i/"^^(Z)^-c- +5/^"M-c" +5/^™^;rc*)) 

(164) 

Being of dimension two, the breaking term Atreafc is soft. As extensively analyzed in |13ll52j . the existence 
of such a breaking does not jeopardize the renormalizability of the theory as well as the introduction of 
meaningful operators which display good analyticity properties. Of course, the fact that the breaking is 
soft is of pivotal importance here. As underlined in [131 153j. the existence of such a breaking is linked to 
the presence of a boundary in field space, namely the Gribov horizon 50. 



^break — I d X 



28 



6.2 Embedding the theory into a more general one 



The standard way of dealing with soft breaking terms is that of introducing them into the starting action 
as composite operators coupled to a suitable set of external sources. In doing so, a more general action 
displaying exact BRST invariance is obtained. Furthermore, the starting action 5*0 and its breaking term 
^break are recovered by demanding that the external sources attain a particular value, usually referred 
as to the physical value (TjElllS]. In order to construct such a generalized action, we introduce external 
sources {V^^ ,U"^,V^l^ ,tj'^l) transforming as 



su: 



ab 

ab 
/if 



jjab 
yab 



sU: 



sv; 



ab 



0, 
0, 



and the invariant action S^^^^^, given by 



(165) 



S- 



Local 



■\ab 9 rabd Ad\ be , j/aci 7-1 
9/ ^jVu + V [^t^ 



ab Ifabdj^d^-bc 



ac 



r^ab 9 rabd /id \ -be 



I MO,cl r^ab 9 rabd /\d\,„be Trac 



jjab _ ^fabd^d^^be ^ gfab^^g^^^ + gf^^A'^^c')^ 



,bc 



abi I 



pdei Ad, 



ab 9 rabd Ad\ - be 



jaoa j^Q 



(166) 



The expression S'Locai is thus recovered from the generalized one •S'L'Jj^g^j when the sources attain their 
physical values: 



V. 



ab 



phys 



-v: 



ab 



phys 



N{N - 1) 



u: 



ab 



fj,u 



phys 



u, 



phys 



0, 



so that 



oinv 



Local 



phys 



'S'Local • 



We end up thus with a generalized invariant action E| 

S'inv = 'S'yM + 5'mAG + ^i^o^al ! 

which reduces to Sq in the physical limit 



phys 



So 



(167) 
(168) 

(169) 
(170) 

(171) 



®The complete invariant action must contain a generalized version of the term Sa in (|17p in order to establish its 
renormalizability. We refer the reader to virhere this analysis is done in detail for the SU{2) case. However, in the 
present section, we are ultimately interested only in the computation of the tree level propagators of the ma^cimal Abelian 
gauge; therefore, we can set q = at this level. 
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The BRST invariance enjoyed by Si^w can be translated in a powerful functional identity, the Slavnov- 
Taylor identity, which is the starting point for the analysis of the all order renormalizability, which we 
shall present separately in a forthcoming work. To derive such identity, we need to introduce more 
external sources coupled to the nonlinear BRST transformations of the fields [54] • To that purpose, we 
notice that the BRST transformation of the off-diagonal component of the gauge field A'^{x) can be split 
in three nonlinear parts: 

sA';, = V; + Ql + n-^, (172) 

jya^_jjab^b^ Qa^^_gjabc^y^ nl = -gr^iAlc\ (173) 

These nonlinear terms can be defined separately by introducing the following set of external sources 

s^^ = T--J7», sTl = sni = Q, (174) 



and writing 



Since 



^£1 = j d'x [j]^ + Qi + Tini + {sQD + (sTep] . (175) 

s'^A"^ = sV + sQ° + s7^'' = , (176) 
one immediately verifies that sS^^^ = 0. The remaining nonlinear transformations can be defined by 

= j d'x [ni, (.4) + L'^ {sc^) + P (sd)] , (177) 

with 

= sV = sL'^ = . (178) 

Therefore, the complete action Sq defined by 

So = 5inv + 5S + 5S, (179) 

obeys the following Slavnov- Taylor identity: 

_ /■,4 U SY.Q 6T,o 6T,o \ 6T,Q dT^p dT^p 6T.oST,o ST,q5T.o 
^ °' ~ J ^ SK^ 5Tl) 5AI 59^^ 5L^ 5c'' 5U be 



0. (180) 



6^- ^ ^ ^'M- 

Besides the Slavnov- Taylor identity, eq. (|180p . the action Sq displays a rather rich symmetry content. 
There exist in fact several symmetries involving the exchange between the Faddeev-Popov ghosts (c", (f ) 
and the localizing auxiliary fields {Lp,(f,uj,Lj), namely 

5^So = 0, 5^So = 0, (i;iSo = 0, (i^So = 0, (181) 
where the operators (5JJ, ^Jl, d^, d^) act on the fields and sources as: 
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• The ^"-transformation: 



• The (^"-transformation: 



• The d" -transformation: 



• The d^-transformation: 



rc -ab X xbc jx 

= \V^^, (182) 



= \u;i, (183) 



di^l = (184) 



JC -a -ac I ^.fabd-bcjl i „fa.bi-bcA 

= -^u +9J ^yC +9J C 

£ 



^ci^a n^fO,bd-bc„c , A^fibi-bcA ■ „2 fabd fdei -be „e 

2 2 
jaMydei-bc^d^e _ ■9_ j^abd fdef -bc^e^f ^ 

^V^''IJ, ~ ''flu 1 

"'i/^ix ~ 2 M"^ ' 



2 ' 
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We can also show that: 

{5l,s} = dl, [5l,s]=dl, K,s]=0, = 0. (186) 

The symmetries of Sq will help us to determine a suitable set of dimension two operators which can be 
consistently introduced in the theory. This will be the subject of the following section. 

6.3 Dimension two condensates 

In this section we shall spend a few words about the subject of the dimension two condensates, which are 
the result of the condensation of local dimension two operators. One should notice that the introduction 
of the horizon function in its localized form, expression (jl30p . entails the introduction of a dimension two 
condensate. In fact, the gap equation (jl06p . implies that the dimension two operator {B^{tp'^ — (p'^)) 
acquires a nonvanishing expectation value, i.e. {B'^^{ip'^ — (p'J^)) ^ 0. An analogous condensate is found 
in the Landau gauge [71 [H [lOl [12l [13] , where the gap equation for the Gribov parameter 7 implies that 
l^^ABC ^^AB _ 0^-8)^1^^ ^ where {4>^^ , <^^^) are the auxiliary fields needed for the localization of the 
horizon function in the Landau gauge. 

Furthermore, in complete analogy with the case of the Landau gauge [101 [T2[ [T3] . other dimension two 
condensates have to be taken into account in the maximal Abelian gauge, see also [l9] for the particular 
case of SU{2). More precisely, the following dimension two operators can be introduced in a way which 
preserves renormalizability of the theory as well as its symmetry content: 

Oa2=A-A-, (187) 

% = <??<'-«-c'^cV (189) 

The operator (jl87p is related to the dynamical mass generation for the off-diagonal gluons, a feature 
which supports the Abelian dominance hypothesis. Its condensation has been established in [2], where 
a dynamical off-diagonal gluon mass has been reported. The ghost operator (jl88p is needed in order to 
account for the dynamical breaking of the SL{2,M) symmetry present in the ghost sector of the maximal 
Abelian gauge: 

5So = 0, 5c^ = c\ 5W = gP^'c^c^ + ^r^^c^c'' . (190) 

Its condensation has been analyzed recently in [55]. Concerning now the third operator, eq. (1189p . we 
notice that it depends on the auxiliary fields ((/?, w, a)). It is in fact needed to account for the nontriv- 
ial dynamics developed by those fields. Besides, one can see that this operator is invariant under the 
transformations (|182p - (|185p . An analogous operator has been found in the Landau gauge [TOl [T2| [T3]. 
where it has allowed to reconcile the Gribov- Zwanziger framework with the most recent lattice data on 
the gluon and ghost propagators [20l [21] . 

Let us briefly show how these operators can be introduced in the theory, by taking the example of 
the operator 0^2 . The other operators can be handled by following an analogous path. Let us introduce 
a BRST doublet of external sources: 

sA = J, sJ = 0, (191) 
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and let us define the following BRST invariant term 



.4 / 1 C 



5*^2 = s I d X [ - X 0^2 — - A J 



d^^E J + - A sOa2 (192) 

where C is a dimensionless constant parameter necessary to account for the ultraviolet divergences aff'ect- 
ing the correlation function 

{Oj,.{x)OA2{y)). (193) 



As discussed in |44j . the parameter C, is uniquely determined by the renormalization group equations. 
Expression (jl92p is thus added to the action Sq, giving 

Si = So + Sa2 . (194) 

Keeping the source J{x) and setting all other external sources to their respective physical valueJ^. we 
can introduce the generating functional W[J] according to 



e 



-W[J] ^ j p^g-So-/d4x(i J0^2-| J2) _ ^-^gg^ 



The vacuum expectation value of the operator is then obtained by differentiating with respect to J: 



5W[J\ 



5J 



= -\{Oa2). (196) 
j=o ^ 



In practice, to solve (jl96p results in a difficult task. A shortcut is usually employed, amounting to make 
use of the Hubbard-Stratonovich field ri{x). To introduce this field in the theory, one inserts the unity 

l=Mjvrje-^^^''<^^'^''-^-''y, (197) 

where A/" is a normalization factor, so that expression (jl95p becomes 

w[j] ^ f^^^^-j:,+jd^x^j 



= So + + ^ 0^2 + 1(0^2)2). (198) 

With W[J] written in this way, we can achieve the following relation 

(0^2) = --(??), (199) 
9 

which easily follows from expression (1198P upon differentiation with respect to the source J. The advan- 
tage of having introduced the Hubbard-Stratonovich field t] relies on the fact that the quadratic term 
in eq. (|195p has been replaced by the term S^, eq. (jl98p . Also, from eq. (jl99p one observes that a nonva- 
nishing vacuum expectation value of the Hubbard-Stratonovich field rj gives a nonvanishing condensate 
(0^2). It remains thus to find out whether the field r] acquires a nonvanishing vacuum expectation value, 



^°We notice that the external sources K^, T|i, i^^, i?^, L"'') and A(a;) carry nonvanishing ghost number, so that their 

physical values vanish. 
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a task which can be accomphshed by evaluating the effective potential corresponding to the action S^. 
A detailed account of the analysis of the effective potential can be found in [H], where a nonvanishing 
vacuum expectation value for rj has emerged. Setting 

r]{x) = {r])+fj{x), 

{fl) = 0, (200) 



we get 



= E„ + / d*. ( + J_ + M o,, + J_ O,, + 1 ip,.r I . (201) 



Introducing thus the gluon mass 



m2 = M , (202) 



we can also write 



Qm rj m rj 



E, = Eo+ / d''x{^ + -^ + — OA2 + -^OA2 + -iOA2y], (203) 



from which one sees that the condensation of the operator 0^2 results in the dynamical generation of a 

2 

gluon mass, i.e. ^ A'^A'Jj^. This term will affect the tree level off-diagonal gluon propagator. In much 
the same way, the other operators Ogxc ^ff will affect the propagators of theory even at the tree 
level, as we shall see in the next section. 



6.4 Tree level gluon and ghost propagators 

In this section we will establish the qualitative behavior of the gluon and ghost propagators by taking 
into account the effects of the restriction to the Gribov region and of the condensation of the dimension 
two operators (jl87p - ()189p . encoded in the following dynamical parameters 

(O^.) ~ , (Oi,^) ~ , {Off) ~ . (204) 

Such parameters will appear in the resulting action as 

/2 
d'^x (^^Oa2+ v' O^xc + /^^ C>ff + "interaction terms" ) , (205) 

where S is the complete action containing all condensates. In order to evaluate the propagators, it is 
sufficient to consider only the quadratic terms of S: 



-'quad 



lim I d X 

/3->0 



abc /tc\ / ab -at 



+ V mfr^) + 2 ^ ) - 



(206) 



where we have already integrated out the Lagrange multipliers {h"',U) and where we have taken the 
physical values of the sources {V^V ,U,U). A further integration over the auxiliary localizing fields 
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(ip, If, u, u) gives the following expression in momentum space 

Sq,,d = lim|^Q^;i(fc)P^^(fe;a)4(-A;) + i4(fc)Q;fU^^ 

-c^ik) W^ik) c\-k) - d{k) 5'^k^ c>{-k) \ , (207) 



where 



7e;=^^^(A;) = 5''\k^ + ^?)-gr^'v\ (208) 
The tree level propagators of the theory are thus given by: 

• The off-diagonal gluon propagator: 
The diagonal gluon propagator: 



(/c2+m2)(/c2+/^2) + ___ 



The symmetric ofF-diagonal ghost propagator: 



(c-(fc)c''(-A;)),,^^ = _ S^' ■ (211) 



/v } ) tiA/rn m — 



The antisymmetric ofF-diagonal ghost propagator: 



^ gjabiyi 

(c (A;)c ( — A;))antisymm = —— — " „2,,2 ! (212) 

(A;2 + ^2)2 + 



The diagonal ghost propagator: 



{c^(ky(-k)) = ^5^^. (213) 



With the exception of the diagonal ghost propagator, eq. ()213p . which exhibits a free behavior, we observe 
that all remaining propagators turn out to be suppressed in the infrared, a result which can be seen as 
a direct generalization of what has been found in the case of SU{2) 
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The off-diagonal gluon propagator deserves a little comment. Notice that it can be written as 




k"^ + m' 



1 



(214) 



where pn{^) represents the deviation from the Yukawa type behavior, and is given by 



PN{k) 



N(N-l) 



(215) 



(A;2 + m2) (A;2 + m2)(A;2 + ^2) + 



(Af-2)g274 
N{N-l) 



The factor p]\f{k) vanishes for = 2, so that the off-diagonal gluon propagator behaves exactly like 
the Yukawa propagator, as obtained in [521 112] for the particular SU{2) case. However, for > 2, it 
seems to deviate from a pure Yukawa behavior, a feature which would be worth to investigate by lattice 
simulations in the relevant case of SU{3). 



7 Conclusions 

In this work we have studied a few aspects of the issue of the Gribov copies in SU{N) Euclidean Yang- 
Mills theories in the maximal Abelian gauge. The so-called Gribov region O has been introduced and 
some of its properties have been established. Summarizing, the region is convex, bounded in all off- 
diagonal directions, and unbounded in all diagonals ones. 

The implementation of the restriction of the domain of integration in the functional integral to the 
region Q has been considered. A careful study of the horizon function has been provided. In particular, 
we have shown that use of Gribov's no pole condition allows us to select only one candidate, given by 
expression (jll7p . A local action from the restriction to the region $7, eq. (jl30p . has been constructed and 
its symmetry content established. 

It is worth mentioning that the results obtained in [32\ Wl\ H9] for the particular case of SU{2) can 
be completely recovered by setting N = 2, which provides a very good check. 

A detailed analysis of the propagators of the theory has been performed. The general case of SU{N), 
N > 2, displays differences with respect to the case of SU{2). In particular, as can be observed from 
expression (j209p . the off-diagonal component of the gluon propagator turns out to be affected by the 
restriction to as well as by the condensation of the operators (jl87p - (jl89p . The diagonal gluon propa- 
gator, eq. (j210p . exhibits a Gribov-Stingl type behavior depending on the Gribov parameter 7 and on the 
parameter p corresponding to the condensation of the operator (jl89p . The symmetric off-diagonal ghost 
propagator, eq. ()21ip . turns out to be dependent from the parameters p and v"^ = v'^v^ while, as expected, 
the antisymmetric off-diagonal ghost propagator, eq. ()212|) . turns out to be directly proportional to the 
parameter stemming from the ghost condensate (Cg^ost) 

All propagators are seen to be suppressed in the infrared. Moreover, they are nonvanishing at A; = 0. 
Although numerical studies of the gluon and ghost propagators in the maximal Abelian gauge have been 
performed only in the case of SU{2), these features seem to be in very good agreement with the most 
recent numerical data [41]. From this point of view, it would be rather interesting to perform a numerical 
study of SU (3) in order to check our prevision. 
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A Details on the calculation of the third order off-diagonal ghost cor- 
relation function 

In order to evaluate the off-diagonal diagonal ghost correlation function {c°'{x)c^{y)), let us start with 
the usual Gell-Mann & Low formula theory: 



where the fields appearing in the right hand side of eq. (j216p are free fields and 5int [c,c] stands for the 
interaction term given by: 



{c^{x)c\y)) = (0|cS(x)c^(y)e-^^-[^"O''=«l|0) , 



(216) 



/ 




(217) 
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To evaluate the aforementioned two-point correlation function we expand the term e till the third 
order, so that 

(c"(x)c^(y)) = (0|c"(x)c^(y)(l-5int + ^5?,t-^5fnt + ---)|0) 

= {Or{x)c\ym - {0\cfix)c\y) S,^,\0) + ^ {0\c^{x)c\y) SfjO) 

-^{0\cf{x)c\y)S?jO) + -- - , (218) 

where 

{0\cf^ix)c\y)\0) = S'^'Goix - y) , (219) 

with 

^^^--y^ = I (04 (220) 

We are interested in terms of order three in the gauge field A or, equivalently, of order g^. These terms 
appear in the expression for S?^^ and Sf^^. Therefore, for the third order correlation function, Gf'{x, y; A), 
one writes 

Gtix,y;A) = l(^^{x)c''{y) J d^xid'^x^ [4g^ T'''' T'''''' f'^^'' A'^Hxi)Ai^X2)Ai^X2) 

+2g^f'^lblcifa2C2i2 j62C2J2^ci (3.^)^i2 1^X2) (xz) + 2g^ f'''^''' f^"^^"^ f2d2i2 j^c^ {xi)Al^ {X2)Ai^ (xz) 

-2g^f-lblc,fC2a2i2fd2b2i2AC^^^^)AC2^^,^)Ad2^^,^)ya^(^^^^ 

+ ^(^<f{x)c''{y) J d''xid^X2d^X3 (8g^n''''T'"''n'''''Ai}ixi)Ai^X2)Al^X3) 

+ I2g^f-Ibrn^a2b2i2fa,b,csj^i^^^^^j^i2^^^-^^C3^^^^ 

+g^r''"^'r''"''n'"''A';i{xi)At'{x2)Amx3)) 

Xe\Xi)dl'c''\xi)e\X2)dl'c''\X2)e\X3)dl'c'''{X3)) ■ (221) 
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Performing the Wick contractions we get 

_^2^3^ai6ici^a2C2»2j62C2J2^ci(3.^)^J^2(^2)^i=^(a;2)+2/ri''l"V"''*'"^ 
X [{dl' {e{x)c''' (Xl)>) (C°^ {X2)c\v)){lt' (Xl)c''^ {X2)) 

+{^{x)c''- {X2)){^' {xi)c\y)) d^' (c"^ {X2)c'" (Xl))] 

y d\id\2d*X3 (8g^r''"''r''"''n'"''^i}{xi)Ai^x2)Ai^x3) 

+12g^r^''^'^r^'^'^r^''^'^Ai}{xi)Ai^X2)Al^X3) + 6g^n'''^ 
+g^n'"'T'''''n'"''A''^^xi)Amx2)A'.Hx3)) 

X [(9^1 {C^xy^ (Xl))) {X2)c\y)) (d^-^ r 1 {X^y-^ {X3)))d:^ r ^ (2.3)0"=^ {X2)) 

\xy' (x-i))) (c"-^ (:C3)c'(y)) {d:^ (c"i (:ci)c'^ (^2))) (c"^ {x2)c'- [x-^) 
\x)c'- {X2))) (c"i {xr)c\y)) (df,' {c^^ {x3)c'' (^1))) 9^'^ (c"^ (:c2)c'- (a;3)) 
'(a;)c''^ (0:2))) {<f' {X3)c\y)) (^^^ (cf ^ (a;i)c'-^ (^3))) (cf = (a;2)c''i (a;i)> 
'(a;)c''^ (X3)>) {c"^ (xi)c''(j/)) (a^^ {c"= {x2)c''' (xi))) dZ^ {^^ {x3)c''^ {X2)) 
\x)c'' {X3))) {X2)c\y))(d:' {c"i {xi)c'' (X2))) {c"=* {X3)c'' (xi))] . 



(222) 



Renaming now some dummy indices, it follows that 

Gt{x,y;A) = i| d\id^X2 [^gT"''' ^^^'^ f^^^'^ A^' {x^)Al^X2)Ai^X2) 

+4g3j<Iliin ja2d2C2^b2d2.2^n(-^..^)^C2(^.2)^j^2(^.^) _ 4^^^ 
+ 2g^f'^lbl'^lf'^202^2fb2C2j2J^c^l(^^^^J^^2(^^,^^^J2^ 

-2pV"'""V''''V'''"''^M (a;i)^?(a;2)vl?(a;2)) 

X [(a;^^ r (:e)c''i (:ei)>) (c"^ (:C2)c''(y)) (c^^ (x-i)c''^ (x'a)} 

+ (<fix)c'^X2)){<^''{x^)c>'{y))d;^i^^X2)c''ix,))] 

-i y d''xid*X2d*X3 (48//"^"^*^"''''= 4' (a;i)X=(a;2)4='(x3) 

+24/ri'i*V"''''V""''"'^Mna;i)X^(K2)A:='(x3) + 24/ri''i'\r^*^"V"'*'^ 

_^24g"/"l''l<=\/°2b2»2 ja3b3i3^ci(,^^-|^i^2 (2.2)^^^(2.3) _^ 

_^^2g:^^«l6l<=l/°2b2»2ja3b3C3^ci(^.^)^^2(2;2)^J3(2.3) + 12//-l''l'=l^ 
^gg3^ai6ici^^a,f„C2^^a3i,3C3^ci(^.^)^c,(^.^)^C3(2.3)^ 

X [(a;^((f (a;)c''Ha;i)>)(c''=(x2)c'(y))(a^Mc''n^i)c'^(a;3)>)9SMc'^^(x3)c''=(a;2)>] . (223) 
It turns out to be convenient to introduce the following quantity 



N{N - 1) 



Grix,y;A). 



(224) 
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Using expression (12190 . we get 

Al{xr)Al{x2)AUx;)) 
X \(dl^Go{x - x^)^Go(x2 - y){dl^Ga{x^ ~ x?S)dl^G^,{xz - 12)] 

+2/r''7'''7''*^7^i)^^(^2)Ai(^2)-2»7"''7^''7""M;(:ci)A^.(x2)A^(a;2) 

+/^'^/''*.^*■^M(^l)^7a;2)Ai(x2)+37"''7''7"*^M(a;l)^:(^2)AL(a;2) 

-g" r""^ fl'^' Al{xr)At{x2)Al{x2)) 

X [(a;;iGo(a; - a;i))Go(x2 - y)Go(xi - ^2) + G„(a; - x2)Go(si - y)a^iGo(s2 - xi)] (225) 
Let us now show that almost all the terms in the double integrals vanish, i.e 

= -r'^r''f''A^^{xi)At{^2)Ai{x2) 

= -r'^r''f''A^^{xi)Ai{x2)At{x2) 

= 0. (226) 



r'7'"7'^"^'4(^i)^^(^2)Ai(x2) = -r''r'''^f''^A^^ix^)Al{x2)Ai{x2) 

^ ^ foMfod^ f'i'^ + rb^r'i^f<ib)A^^(x,)Alix2)Aiix2) 

_Jabi Jabd — Q 

= r''r'''r'''A^^{xi)Ai{x2)Aiix2) 

= -r''7""V^4(^l)^^(x2)Ai(x2) 

= 0, (227) 



r''r''f'''''A^^ix,)Alix2)At{x2) 



r'''f'''r''^A';^ixi)Al{x2)Ai{x2) 



-r''nf'"'^A^,{xMii^2)At{x2) 
-r'^f'''r''A',{xMt{x2)Ai{x2) 
-r'^f'''r''A',{xMi{x2)At{x2) 

0, (228) 

-r''r'''f''^Ai{x^)A'^{x2)A^Ax2) 
-r'7"'V*^^(^i)^i(x2)4(^2) 
-r''r'''f'''A'^{xi)A'A^2)A^A^2) 

0, (229) 
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= -r''T'f''''Al{xi)At{x2)Al{x2) 

= . (230) 

The only term that does not vanish is f°''"'f'"^f"''^'A'j^{xi)Al{x2)Ai{x2). Thus, we can finally write Q^^^ 
as 

g'^'\x,y;A) = d\^d''x2d\s(8g'r'^f'^^r''Al,ixr)Ai{x2)AUx3) 

+57"'7''7'"^'^m(^-i)^7^-2)^^(^3)) 

'^ N{N- 1) / '^"^i'^^^^ g^r'"'f'"^T'''Al,{xi)Al{x2)Al{x2) 

X [(^a^iGo(a; - xi))Go(a;2 - y)Go{xi - X2) + Go(x - a;2)Go(xi - y) dl^Go{x2 - xi)] . (231) 
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